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ABSTRACT 

The splitting of photons 7 ^ 77 in the presence of an intense magnetic field has 
recently found astrophysical applications in polar cap models of gamma-ray pulsars 
and in magnetar scenarios for soft gamma repeaters. Numerical computation of the 
polarization-dependent rates of this third order QED process for arbitrary field strengths 
and energies below pair creation threshold is difficult: thus early analyses focused on 
analytic developments and simpler asymptotic forms. The recent astrophysical interest 
spurred the use of the S-matrix approach by Mentzel, Berg and Wunner to determine 
splitting rates. In this paper, we present numerical computations of a full proper-time 
expression for the rate of splitting that was obtained by Stoneham, and is exact up 
to the pair creation threshold. While the numerical results derived here are in accord 
with the earlier asymptotic forms due to Adler, our computed rates still differ by as 
much as factors of 3 from the S-matrix re-evaluation of Wilke and Wunner, reflecting 
the extreme difficulty of generating accurate S-matrix numerics for fields below about 
4.4 X 10^^ Gauss. We find that our proper-time rates appear very accurate, and exceed 
Adler's asymptotic specializations significantly only for photon energies just below pair 
threshold and for supercritical fields, but always by less than a factor of around 2.6. We 
also provide a useful analytic series expansion for the scattering amplitude valid at low 
energies. 
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1. INTRODUCTION 

The exotic quantum electrodynamical (QED) process of the sphtting of a photon into two 
(7 ^ 77) in the presence of a strong magnetic field has recently been of interest in the study 
of gamma-ray pulsars and soft gamma repeaters (SGRs). Many radio pulsars are known to have 
spin-down fields in excess of 10''^^ Gauss, with about a dozen exceeding 10^^ G, including the 7 -ray 
pulsar PSR1509-58. This source is unusual among the seven known 7 -ray pulsars in that it has 
a hard power-law spectrum that turns over sharply around 1 MeV (Ulmer et al. 1993; Hcrmsen 
et al. 1997), whereas other pulsars (e.g. the Crab; see Nolan et al. 1993) whose fields are nearer 
10^^ G have emission extending beyond 1 GeV. A question then naturally arises: does the high 
field of PSR1509-58 cause its spectrum to differ so much from that of the other 7 -ray pulsars? It 
appears that magnetic photon splitting, in concert with the more familiar process of single-photon 
pair creation, may provide the answer (Harding, Baring & Gonthier 1996, 1997) in the context of 
polar cap models, inhibiting emission above a few MeV for a reasonable range of cap sizes. Photon 
splitting may also be quite important in the most recently discovered (Ramanamurthy et al. 1996) 
gamma-ray pulsar, PSR0656-I-14, which has a spin-down field of 9.3 x 10^^ Gauss, between that 
of the Crab and PSR1509-58. 

Also of interest to high energy astrophysicists are soft gamma repeaters (SGRs), the transient 
sources that are observed to have sporadic periods of high 7 -ray activity. There are three known 
repeating sources, all with sub-second durations and soft spectra ( kT ~ 30 keV: for SGR1900-I-14, 
see Kouvelioutou et al. 1993). The identification of radio (SGR1806-20: Kulkarni and Frail, 1993) 
and X-ray (e.g. Mar 5, 1979 repeater: Rothschild, et al. 1994) counterparts to SGRs spawned 
much interest in these objects. It has been suggested that SGRs are magnetars, neutron stars with 
extremely high fields, in excess of the quantum critical field strength Bcr = mlc^ /efi = 4.413 x 10-^^ 
Gauss (when the cyclotron and electron rest mass energies are comparable); the spin-down estimate 
for the Mar. 5, 1979 repeater is i? ~ 6 X 10^^ Gauss (Duncan and Thompson, 1992), assuming 
its age is that of N49, its associated supernova remnant. Such fields far exceed those found in 
most radio pulsars, and can perhaps be generated by enhancement due to dynamo action (Duncan 
& Thompson, 1992). Baring (1995) proposed the potential importance of photon splitting for 
softening emission spectra in SGRs via a splitting cascade to reproduce the observed emission 
spectra, yielding an estimate of B ^ 2 x 10^^ Gauss for SGR1806-20. The observed stability of 
the spectra from outburst to outburst argues in favour of an equatorial emission region (Baring 
& Harding 1995; Harding & Baring 1996). Splitting can also be quite infiuential in SGR spectral 
models in the absence of cascade formation (e.g. Thompson and Duncan 1995). 

To date, astrophysical models that invoke splitting in the environments of gamma-ray pulsars 
(Harding, Baring &: Gonthier 1996, 1997; Chang, Chen and Ho 1996) and SGRs (e.g. Baring 
1995; Baring & Harding 1995; Thompson and Duncan 1995; Harding & Baring 1996; Chang et 
al. 1996) use relatively simple approximations to the splitting rates that were derived by Adler 
et al. (1970, see also Adler 1971; Bialynicka-Birula &; Bialynicki-Birula 1970), using effective 
Lagrangian techniques, which are valid for low photon energies emeC? and fields satisfying eB ^ 
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Bcr ■ Restriction to this asymptotic regime has until recently been necessitated by the total lack 
of more general evaluations of splitting rates that are sufficiently amenable for use in astrophysical 
models: formal expressions for the third-order QED process of 7 — >■ 77 are very complicated. This 
has been a real deficiency, since the eB ~ Bcr regime can easily be realized in both SGRs and 
pulsars. This paper addresses this limitation by presenting numerical determinations of splitting 
rates for arbitrary field strengths and photon energies below the pair creation threshold of 2me(? , 
as obtained from the general proper-time QED formulation (a technique due to Schwinger) of 
splitting presented by Stoneham (1979). 

In this paper, our numerical computation of Stoneham's splitting rates focuses on the polariza- 
tion mode -L^llll that is permitted in the limit of weak dispersion by energy-momentum kinematic 
selection rules (see Adler 1971); ±^_L_L and ||^-L|| are the other two modes permitted by QED in 
the non-dispcrsivc limit, and results for them are similar to those presented here (Baring, Harding 
and Weise, in preparation, hereafter BHW97). Here _L and || denote the orientation (perpendic- 
ular and parallel to, respectively) of the photon electric vector relative to the plane containing its 
momentum and the field line. Our computed rates vary continuously from the low energy limit, 
and reproduce widely-accepted analytic asymptotic forms at both low and highly supercritical 
field strengths. The numerical computations presented here therefore appear extremely reliable for 
use in astrophysical applications, and concur with the recent numerical work of Baier, Mil'shtein, 
Sz Shaisultanov (1996), which resulted from their alternative Schwinger-type analysis of photon 
splitting. Note, however, that our rates differ dramatically from the numerical evaluation of the 
S-matrix determination of the splitting rate by Mentzel, Berg and Wunner (1994, sec also Wunner, 
Sang and Berg 1995) that is orders of magnitude greater than the earlier proper-time determina- 
tions, and differ by as much as a factor of around 3 from the recently-updated S-matrix numerics 
of Wilke and Wunner (1996; see also the new Erratum to Wunner, Sang and Berg 1995), which still 
show an inability to produce accurate results at subcritical field strengths, the domain applicable 
to pulsar models. This discrepancy underlies the inherent difficulty in computing S-matrix rates 
for relatively low fields, and contrasts the reliability of the proper-time computations presented in 
this paper and also by Baier et al. (1996). 



2. PHOTON SPLITTING RATES AND SPECTRA 

Magnetic splitting 7 77 is a relatively recent prediction of QED. A decade of controversy 
followed the earliest calculations before Bialynicka-Birula & Bialynicki-Birula (1970), Adler et 
al. (1970) and Adler (1971) performed, via an effective Langrangian technique, the first correct 
evaluations of its rate, including asymptotic forms in the limit of photon energies well below pair 
creation threshold, which varied as B^ when B ^ Bcr ■ Their rate determinations neglected 
dispersion in the birefringent, magnetized vacuum, though the effects of such dispersion on the 
possibility of energy and momentum non-conservation and associated polarization selection rules 
was discussed in detail by Adler (1971). The early controversy was fueled by the inherent difficulties 
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in calculating the rates of this third order process by standard QED techniques, and the physics of 
photon splitting at high magnetic field strengths is still the subject of debate. 

The most general calculation of photon splitting rates was performed by Stoneham (1979), who 
used Schwinger's proper-time method to determine formal expressions for the splitting rate for e 
below pair creation threshold and arbitrary field strengths. This analysis fully included dispersive 
effects due to magnetized vacuum in formal expressions for the rates, leaving the rates for the 
different polarization modes of splitting as unmanageable nine-dimensional integrals. Fortunately, 
in the limit of zero vacuum dispersion, Stoneham's formulae specialize to amenable triple integrals 
that involve only moderately complicated combinations of elementary functions (see Adler 1971, 
Papanyan & Ritus 1972, and Baier, Mil'shtem, & Shaisultanov 1986 for alternative presentations, 
and the recent compact results of Baier, Mil'shtem, & Shaisultanov 1996, hereafter BMS96; Adler 
and Schubert 1996). We have analytically developed these triple integral expressions (see BHW97) 
to facilitate our numerical purposes, eliminating leading-order terms (linear in photon energies) that 
cancel exactly, leaving scattering amplitudes that are cubic in photon energies. Such cancellations 
plague numerical evaluations of splitting rates, and can be a principal source of error. Hence their 
elimination via an algebraic development is highly desirable. Readers interested in the details of our 
analytic reduction of Stoneham's (1979) expressions can refer to the more extensive presentation 
in BHW97, where the full expressions that are computed in this paper are exhibited. 

We note that while such analytic developments are expedient and very useful, it is possible to 
proceed by numerically subtracting off cancellations in a well-chosen manner. This is the approach 
discussed in Adler and Schubert [1996; see their Eq. (17)], who have developed codes to compute 
the integrands (i.e. double-integrals) of the triple integrals composing the scattering amplitude. 
Their codes, which are publicly available on the World Wide Web (the address is listed in Adler 
and Schubert 1996), demonstrate numerically the equivalence of the integrands that emerge from 
the different formalisms of Adler (1971), Stoneham (1979), BMS96 and also Adler and Schubert's 
presentation. Unlike the results presented in this paper and the work of BMS96, Adler and Schubert 
(1996) do not compute the full triple integral for the amplitude, a computationally demanding 
extra step that is necessary to obtain results of use for astrophysical applications. The numerical 
algorithm used to obtain the results of this paper involves a mixture of Simpson's rule and Gauss- 
Laguerre quadrature to obtain accurate differential and total splitting rates. 

For quite general photon energies and magnetic field strengths, the differential attenuation 
coefficient (i.e. the rate divided by c) of the photon splitting mode -L— ^|||| can be written in the 
form (BHW97) 



for an initial photon of energy e and produced photons of energies e\ and £2 = e — £1 . Here 
M. denotes the scattering amplitude, scaled according to Stoneham's notation; other polarization 

modes assume a similar form. Hereafter, dimensionless units will be used, with B being expressed 
in units of i?cr and photon energies e (initial) and £\ and £ — £1 (final) being scaled by mgC^ . 




(1) 
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In equation (||), ai is the fine structure constant, A = 3.86 x 10~^^ cm is the electron Compton 
wavelength over 2tt , and 9 is the angle the photon momenta make with the field lines; equation (|T|) 
is strictly valid only in the zero-dispersion limit, for which the initial and final photon momenta are 
collinear. The triple integral expression we obtained for the scattering amplitude A4 , somewhat 
lengthy for presentation here, has an integrand consisting purely of terms that are cubic combina- 
tions of e — El and ei , multiplied by an exponential that is quadratic in these energies, all divided 
by the common factor eei{e — £i) ■ 

When the incident photon energy is low, namely e <C 1 (i.e. 511 keV), these energy depen- 
dences in the integrand cancel exactly, the exponential loses its energy dependence and dramatic 
analytic simplification is possible. The scattering amplitude then becomes dependent only on the 
field strength: 



3 + 2s^ , s cosh s 



for arbitrary B . The form of A4i{B) is just that given in Adler (1971) and Eq. 41 of Stoneham 
(1979), and is that used in the recent astrophysical models of SGRs (e.g. Baring and Harding 1995) 
and gamma-ray pulsars (e.g. Harding, Baring and Gonthier 1996, 1997). In the limit of B <^ 1 , 
Mi{B) « 26/315 , while in the limit of S > 1 , Mi{B) ^ 1/{6B^) . For arbitrary B and energies 
well below pair creation threshold, a series expansion of the scaled scattering amplitude 
in terms of photon energy was obtained by BHW97: 



= Mi(B) + 



(^el + el) Mu{B) + eie2Mi2{B) 



sin^ 6* + . . . , (3) 



where the produced photons assume dimensionless energies ei and 82 = £ — £i , and 

„2 - 



3 f ^ I ^ \ coshs I / 19 s \ coshs 1 

SsinhS ~ I 47 30 j ihSI + I 64i + 24 i ( 



and 



15 / 13 s \ cosh s I 1 5s \ cosh g I 

^32sinh4s " \8? ~ 30 j ihS7 + l^TBi " 24 j j ' ^ ' 

The expansion in Eq. (|3|) is a result that is quite useful to pulsar modellers that can be reliably 
used for e ^ 0.5 . For B ^ 0.3 , this series reduces to 

A^^^llll |^ + ^^^(205e? + 19eie2 + 205e^) sin^^ , ei, £2 < 2/sin^; S ^ 0.3 . (6) 
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Recently, Baier, Mil'shtem, & Shaisultanov (1996) have derived the high field limit (i.e. B ^ 1) 
for splitting at arbitrary energies below pair creation threshold (e < 2); in the notation used here, 
their result becomes (for 9 = 90° ) 

A^^^llll « 75T^ I arcsinff ) + arcsin(^) - e] , e < 2 , (7) 

for £2 = e — ei . We have reproduced this formula identically from Stoneham's (1979) formulae 
and also our algebraic reduction of his expressions (BHW97). For photon energies far below pair 
creation threshold, equation (0) approaches 1/(65^), the high field limit of Aii{B) , and the 
series expansion in energy maps over to the high field limit of Eq. (^), for which A^n ~ —Ai\2 ~ 
1/(30-6^) . The inverse cubic dependence of this result on the field implies an attenuation coefficient 
that is virtually independent of B in highly supercritical regimes. 

Total splitting attenuation coefficients are obtained by integrating equation (Q) over produced 
energies < ei < e. For e <C 1 and S <C 1 , A^_l_+|||| A^i(O) = 26/315 and this integration is 
trivial, leading to a total attenuation coefficient of 

= I'T^^Wwdei ~ To = (^)^ B^ sin"^ 9 « 0.116 e^sSin^ ^ cm-^ (8) 

INI 7o 607r A Voi<3/ 

In the limit of highly supercritical fields ( S ^ 1 ), no analytic evaluation of the total splitting rate 
[oc efel times the square of the formula in equation (^)] is apparent, as extensive development 
leaves some intractable terms. However, the approximation 

" 3?^X--'nWT-& + St-m8 + W)' 0<.. = esm^<2; B » I. (9) 

is accurate to better than 1.5% at all energies below pair creation threshold, and suffices for nu- 
merical purposes. The asymptotic formulae represented by equations (^) and (|^ and also the 
approximate result in equation (|9|) provide strong checks on the accuracy of the numerical compu- 
tations presented here. 

The computed attenuation coefficients for the splitting mode -L— >|||| obtained by integrat- 
ing equation ([l|) over produced photon energies ei are displayed in Figure la, for different field 
strengths, together with the magnetic pair creation rates for the _L polarization (discussed in 
Daugherty and Harding 1983; Harding, Baring and Gonthier 1997). Two characteristic properties 
of QED radiation processes in strong fields are immediately apparent. First, photon splitting, like 
pair creation, is a strongly increasing function of both the magnetic field strength (when B ^ 1) 
and also of the energy of the incident (absorbed) photon. Second, at energies near pair threshold, 
as a third-order QED process, photon splitting has attenuation coefficients several orders of magni- 
tude smaller than pair creation, a first-order process. Another striking feature of Figure la is how 
closely the coefficients resemble a power-law, even for highly supercritical fields; equation (|9|) 
clearly reveals that the high B asymptotic limit (the dashed curve in Fig. la) has at most modest 
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(i.e. less than a factor of around 2.6) deviations above the dependence near pair threshold. 
Further, Adler's (1971) simplest analytic form (oc e^B^) is quite accurate below pair creation 
threshold for B ^ 0.2 . Deviations from behaviour are illustrated in Fig. lb, where the ratio 
Tj^^lll/To is plotted as a function of incident photon energy. For the different field strengths, 
quasi-horizontal portions of the curves indicate dependence; clearly from the figure, this simple 
asymptotic form is obeyed over a wide range of phase space. Fig. lb also demonstrates the impres- 
sive agreement between the numerical results obtained here and those of BMS96 and Adler (1971). 
The intersections of the various curves with the ordinate axis roughly defines corresponding values 
of (315/26)Ali(5). 

Figure 2 exhibits differential photon production spectra, i.e. A^5_-+|||| ' normalized to unity. 
A remarkable property emerges: except for very near pair threshold, or for highly supercritical 
fields, the spectra closely approximate the low energy asymptotic limit of ef (e — ei)^/30 . Magnetic 
broadening of the spectrum (i.e. with increasing B), & feature of other strong field QED radiation 
processes, does however become very significant near pair creation threshold. These numerical 
results are also in accord with the recent alternative presentation of BMS96. 

3. DISCUSSION 

Motivation for this numerical computation of proper-time results has been enhanced by a 

new result on the rates of photon splitting: Mentzel, Berg & Wunner (1994, hereafter MBW94) 
presented an S-matrix calculation of the polarized splitting rates. While their formal development 
is equivalent (Weise, Baring and Melrose 1997) to an earlier S-matrix formulation of splitting in 
Melrose & Parle (1983a,b), their presentation of numerical results (see also Wunner, Sang & Berg 
1995, where the implications of their calculations were discussed in astrophysical contexts) was in 
violent conflict with the splitting results obtained by Adler et al. (1970), Adler (1971), Bialynicka- 
Birula & Bialynicki-Birula (1970), Papanyan and Ritus (1972), Stoncham (1979) and a number of 
other more recent works. The numerical rates of Mentzel, Berg & Wunner (1994) for B = 0.1 
and B = 1 are larger by orders of magnitude than our computations (and also those of BMS96). 
Their S-matrix rates exhibit a weak energy dependence well below pair threshold, and also exceed 
the pair creation rate at threshold ( e = 2 ) for S = 0.1 . Wilke and Wunner (1996) have very 
recently retracted the earlier S-matrix rates (sec also the recent Erratum for Wunner, Sang and 
Berg 1995), confirming a coding error in the numerics of MBW94. Their numerical re-evaluation, a 
sample of which is depicted in Figure la, is much closer to the proper-time results, but still differs 
significantly (by as much as a factor of 3.3 at e = 0.1 , B = 1, and 45% at e = 1.9, B = 1 
for the illustrated cases) from our determinations, unless the field is highly supercritical. The S- 
matrix and proper-time techniques should produce equivalent numerical results, and indeed have 
done so demonstrably in the case of pair production (see Daugherty & Harding 1983; Tsai & Erber 
1974). The proper-time analysis, though non-trivial, is computationally much more amenable than 
the S-matrix approach, and has been reproduced in the limit of -C 1 by numerous authors. 
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Clearly, the work of Wilke and Wunner (1996), which does not exhibit any results below about 
0.6i?cr , still demonstrates an inability to produce suitably accurate results at critical and sub- 
critical field strengths. This point is not made clear in the recent Erratum to Wunner, Sang 
and Berg (1995). These inaccuracies can dramatically impact astrophysical applications to pulsar 
models (e.g. Harding, Baring and Gonthier 1997), where factors of 3 become very significant. Such 
discrepancies underline the inherent difficulty in computing S-matrix rates for relatively low fields, 
where many Landau level quantum numbers for the three intermediate states must be summed 
over (e.g. sec Weise, Baring and Melrose 1997), and contrasts the reliability of the proper-time 
computations presented in this paper, and by BMS96, that emanate from relatively simple integrals. 

In conclusion, this paper presents the major features of our numerical computations of proper- 
time expressions for photon splitting attenuation coefficients and spectra. Our numerical results 
are quite consistent with several analytic limits obtained by various groups, including the very 
recent approach of Adler and Schubert (1996), comparable to the recent numerical computations 
of BMS96, and well-behaved in the light of physical intuition related to strong field QED processes; 
they still differ significantly from and represent a very useful improvement over the S-matrix re- 
evaluation of Wilke and Wunner (1996), which becomes problematic at field strengths below ~ 
4 X 10^^ Gauss. We believe that our proper-time computations can be reliably and efficiently used 
in astrophysical models. 

We thank Stephen Adler and George Pavlov for reading the paper and for comments helpful 
to the improvement of the manuscript. This work was supported through Compton Gamma-Ray 
Observatory Guest Investigator Phase 5 and NASA Astrophysics Theory Program grants. MGB 
thanks the Institute for Theoretical Physics at the University of California, Santa Barbara for 
support (under NSF grant PHY94-07194) during part of the period in which work for this paper 
was completed. We thank Ramin Sina for making his pair creation code available to us for use in 
the preparation of Figure la. 
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Fig. 1. — (a) The computed photon splitting attenuation coefficients for polarization mode -L— >|| || , 
as functions of the incident photon energy e, for field strengths B = 0.1, 1, 10 (solid lines below 
2me(? , in ascending order), and the asymptotic rate resulting from the insertion of Eq. in 
Eq. (HI), labelled as -B ^ 1 (dashed curve). Photons are assumed to propagate orthogonally to the 
field lines {9 = 90° ). For comparison, the pair creation rates (Daugherty and Harding, 1983) for 
_L photons at two different field strengths are depicted above 2me(? ■ The filled circles denote the 
S-matrix determination of Wilke and Wunner (1996) for B = 1 . 
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Fig. 1. — (b) The ratio of the total attenuation coefficient to the low-energy, low-magnetic 

field asymptotic limit Tq [see Eq. (^] as a function of energy for different field strengths B , as 
labelled. The open squares represent ratios determined from the computations of Adler (1971), 
with the two points near e = 0.1 being the e = results from Adler's Fig. 8. The filled triangles 
are derived from the numerical evaluations in Fig. 2 of BMS96. Impressive agreement between our 
results and those of Adler (1971) and BMS96 is evident. 
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Fig. 2. — The computed photon spUtting differential rates, normahzed to unity and expressed as 
functions of the energy ei of one of the produced photons, for polarization mode -L— >|||| . The top 
three curves are for field strengths B = 0.1 , B = 100 and B = 10^ for initial photon energies 
e = 1 (i.e. 511 keV), with the broad distribution corresponding to e = 2 and B = 10^ . The 
B = 0.1 curves closely resemble the asymptotic low energy form of ef(e — £i)'^/30, which can be 
deduced from Eq. (||), while the B = 10^ distributions closely resemble the shape of the square of 
the high field limit in Eq. (|^. 



